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We derive and evaluate expressions for the low temperature dc equilibrium tun- 
neling conductance between parallel two-dimensional electron systems. Our theory is 
based on a linear-response formalism and on impurity-averaged perturbation theory. 
The disorder broadening of features in the dependence of tunneling conductance on 
sheet densities and in-plane magnetic field strengths is influenced both by the finite 
lifetime of electrons within the wells and by non-momentum-conserving tunneling 
events. Disorder vertex corrections are important only for weak in-plane magnetic 
fields and strong interwell impurity-potential correlations. We comment on the basis 
of our results on the possibility of using tunneling measurements to determine the 
lifetime of electrons in the quantum wells. 
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I. INTRODUCTION 



Tunneling of electrons has long been studied both as a fundamental manifestation of 
the quantum nature of matter and as a powerful probe of many different electronic systems 
It has only recently become possible to measure the tunneling conductance between 
two-dimensional electron systems (2DES's) in semiconductors. The 2DES's occur in 

GaAs quantum wells which are separated by thin AlAs tunneling barriers The strong 
constraints imposed by energy and momentum conservation in 2D can potentially lead to 
very sharp features in the dependence of the equilibrium tunneling conductance on the sheet 
densities in the quantum wells and on the strength of an in-plane magnetic field. As we dis- 
cuss below, tunneling conductance studies probe the properties of 2DES's in a way which is 
qualitatively different from conventional in-plane transport studies. In this paper we present 
a theory which provides a framework for the detailed interpretation of these experiments and 
apply it to study the influence of disorder on the tunneling conductance. We find that these 
experiments can provide a direct measure of the quantum lifetime of electrons in the 2DES's. 
When combined with measurements of the 2DES transport lifetime from in-plane transport 
experiments this provides information on the nature of the disorder scatterers in the sys- 
tem. We also show that vertex corrections to the naive result for the impurity-broadening 
of features in the tunneling conductance are important only for strong correlations between 
the disorder potentials in the separated 2DES's and for weak in-plane magnetic fields. In 
Section II we briefly review the linear response formalism which is the basis of our theory. 
In Section III we present results for the dependence of tunneling conductance on the density 
difference of the 2DES's and the strength of an in-plane magnetic field in the absence of 
correlation between the disorder potentials in the two layers. In the ideal disorder-free limit 
the tunnelling conductance can be evaluated analytically and diverges at certain magnetic 
field strengths. We examine how these divergences are suppressed by disorder and discuss 
the possibility of using the resulting tunneling conductance curves to measure the quantum 
lifetime of electrons in the double quantum well systems. In Section IV we examine vertex 
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corrections to the naive expression for the tunnehng conductance in disordered systems. We 
find that when the disorder in the two 2DES's is strongly correlated features in the conduc- 
tance curves may sharpen. Vertex corrections are important only at weak in-plane magnetic 
fields. Finally in Section V we briefiy summarize our results. 

II. LINEAR RESPONSE FORMALISM 

We restrict our attention here to the case of zero perpendicular magnetic field so that 
states of the 2DES's may be described in a plane-wave-state basis. Our interest here is in 
studying the dependence of the tunneling conductance on experimentally tunable parameters 
like the 2DES layer densities and the strength of an in-plane magnetic field rather than on its 
absolute value. We therefore describe the tunneling phenomenologically using a tunneling 
Hamiltonian 

Here we use / and r to label states localized in left and right 2DES's respectively. The 
remaining terms in the Hamiltonian may include the interaction of electrons in either 2DES 
with a disorder potential and interactions between electrons in the same or in different 
2DES's but are assumed to commute separately with the total number operator for each 
2DES, Ni^r- Following a familiar line we note that the operator describing the net current 
flowing from the left 2DES to the right 2DES is given by 

J = -ej^Nr = ejNi = te[Nr,HT]/fi 
= f E ih',kCl,^^C^,,., - h.c.) (2) 

fc'fc,(7 

Since the expectation value of / is zero in any eigenstate of the Hamiltonian it follows that 
no current flows in equilibrium. We evaluate the current which flows in response to a electric 
potential difference V{t) = V exp{iujt) between the two wells. The perturbing term in the 
Hamiltonian is 
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Hr{t) = -eV{t){Ni-N,)/2. (3) 

Using finite temperature linear response theory we find tliat in tlie limit ^ the tunneling 
conductance, G = (/) /V is given by 

G = hm '^^l^tlinl (4) 

a;^0 hu; 

where K{z) is the analytic continuation of the Matsubara current-current correlation func- 
tion: 

Kizun) = - r dr{Tr[iir)m]). (5) 
Jo 

Eq. (H) is the basic expression we use to evaluate the tunneling conductance. 

III. UNCORRELATED ELECTRON LAYERS 

We first apply Eq. (j^) to the case where the two 2DES's can be considered to be decoupled 
in the absence of the perturbation p|. It then follows from the translation invariance of each 
well and the definition of the Matsubara Greens function |^ that 

^ k,k' «'^» 

+ gj:^^{iujn)gj^,^iiiujn - i^)] (6) 

(We have assumed spin degeneracy.) The frequency sum may be evaluated by contour 
integration methods with the result that 

Here Aj:^{ijj) = — 2ImG^j^(u; + irj) is the spectral weight of the exact finite temperature 
Greens function in layer m. In the absence of interactions and disorder Aj^^ = 2Tih5{e — 
em{k)) and Eq. (|^ reduces to the Fermi golden rule expression for the tunneling conductance. 
Deviations from this limit probe many-body effects in the 2DES's. We see from Eq.( |^ that 
when correlation between the two-layers may be ignored the tunneling experiments measure 



the one-body Greens function of the individual layers, i.e. they measure the spectrum of 
states created when an extra electron is added to the ground state of a system. In contrast 
in-plane transport experiments measure a two-body Greens functions of the layers, which 
reflect the correlations of particle-hole excitations of a system. In many-body systems the 
information available in these two experiments is complementary. 

The numerical results we present below are all for the ideal case where the in-plane com- 
ponent of the momentum is conserved during tunneling, i.e. where ^ = 15^, ^. Physically 
this limit corresponds to the case where the tunneling amplitude is uniform over the area of 
the 2DES's and the comparison of our results with experiment |^ discussed below indicates 
that this uniformity is achieved to a remarkable degree. The results presented below also 
do not account for electron-electron interactions which do not appear to play any important 
role |jlO[ experimentally. We will restrict our attention to T = 0. 



When momentum is conserved during tunneling it follows from Eq. ( ^ that in the absence 
of interactions and disorder G = unless the 2DES's have identical densities. For this reason 
it is useful to study the dependence of G on both the layer densities of the 2DES's and on 
the strength of an in-plane magnetic field, B = Bey. If we neglect the finite thickness of 
the electron layers, the effect of B on the electronic wavefunctions is just to add a gauge 
transformation phase factor which shifts the origin of momentum space in the two layers. 
For example if we choose the gauge in which the vector potential is A = Bze^, it will shift 
the origin of momentum space by —Zm(ix/C-'^ in a layer located in the plane z = z^- (Here 
i = {hc/eBy/'^ is a magnetic length.) With these shifts the single-particle energies in layer 
m are given by 

eUk) = E^^ + ^[{K + zJfy + kl] (8) 

where is the subband energy in layer m. (The layer density difference is related to the 
subband energy difference by Ur — rii = {Ei — Er)/i'o where uq = m/Trh^ is the 2DES density 
of states.) In Fig.( |^) we show the tunneling conductance versus in plane-field strength for 
the situation studied experimentally in Ref.( 0). The dashed line shows the result without 



any disorder which can be obtained analytically: 

G AeH^m*^ 



(9) 



when both terms in square brackets are positive and is zero otherwise. Here K = {zr — 
-2/)/^^ ~ -B is the amount by which the two Fermi circles are displaced as illustrated by the 
inset of Fig.( |l|), A is the area of the 2DES's and kp^ is the Fermi radius in layer m. Note 
that for equal densities in the two layers [kpi = kp^) G diverges as B^^ for i? ^ but only 
as {Be — B)^^/'^ for B ^ Bc = 'n\zr — Zi\/kp^Q. {Be is the maximum field for which tunneling 
can occur and $o is the magnetic flux quantum.) 

These divergences of G are suppressed by disorder. The solid line in Fig.( |l]) shows 
results obtained using the Born approximation expression for the disorder-averaged Greens 
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G^. Ju ± IT,) = {uj- h-hUk) ± ^r'. (10) 

Here r is the quantum lifetime of the electrons. The results in Fig.( |1]) were calculated with 
Tep/h = 15. This level of disorder is important in suppressing the divergences ai K = 2kF 
and especially a.t K = but has relatively little effect on G at other values of K. For equal 
density 2DES's the maximum conductance reached at zero in-plane magnetic field {K = 0) 
is proportional to r: 

G(A- = 0) AV„ ^^^^,^j^j^2A%r ^^^^ 



Here Ab{u!) = {1/t)/{uj'^ + (2/r)^) is the Born approximation Lorentzian spectral function 
and we have assumed that rep ^ 1. This suggests the possibility of using the shape of the 
tunneling conductance curve to measure the electron lifetime. For example, assuming the 
G{K = kp) is nearly independent of disorder gives 

Eq.( ^) provides a lower bound on r since the divergence of G{K = 0) will also be suppressed 
by finite temperatures and also possibly by non-momentum-conserving tunneling events and 

6 



macroscopic inhomogeneities. Another estimate of rep/ k may be obtained by examining the 
dependence of G on the density difference between the two layers. In the Born approximation 
we find that at K = 

^ ^ h\l + 4{{ni-nr)/{ni + nrmeFT/hy) ^ ^ 

In Fig.( ^) we show results for the layer density dependence of the conductance at K = 
and at several non-zero values of in-plane magnetic field. According to Eq. ( |13]) the relative 
density difference at which G{K = 0) is reduced to half its equal density value is epr/h. 

Comparing with the experimental data of Eisenstein et al. Eq.( |12D gives the estimate 
epr/h ~ 5 while the half- maximum density- difference value gives the estimate epr/h ~ 10 
for their sample. These values should be compared with the measured transport lifetime for 
the sample of Ref . [Q, , Ttrep/fi ~ 250. In high-mobility quantum well structures, scattering 
is dominated by ionized impurities which are spatially separated from the 2DES's. These 
impurities produce small-angle-scattering of electrons because the Coulomb potential they 
produce is smooth in the 2DES layers. Small angle scattering is ineffective in limiting 
the transport lifetime. For disorder dominated by remote ionized donors large values of 
Ttr/T ~ {kpdsY therefore expected and it is widely appreciated that the transport 



lifetime no longer effectively characterizes the degree of disorder in the sample, [dg is the 
distance between the 2DES and the ionized impurities.) The estimates of rep/h derived 
from the tunneling date above are in reasonable accord with expectations, although the value 
obtained from Eq.( |l^) is probably an underestimate. The peak value of the conductance is 



likely to be reduced by large length scale inhomogeneities in the samples |]T3[ which we are not 
adequately described in the Born approximation. Any momentum non-conservation during 
tunneling would also lead to further broadening of the features in these curves. The fact that 
the rep/h values inferred without accounting for momentum non-conserving tunneling events 
are so close to expected values demonstrates that momentum non- conserving tunneling 
events are remarkably infrequent in the experiments of Ref.( [Q). This property allows the 
tunneling experiments to probe the properties of the 2DES's in more detail than would 



otherwise be possible. 



IV. VERTEX CORRECTIONS 

Now we will turn to the discussion of the case where correlations exist between the 
disorder potentials in the two 2DES's. It is useful to begin with a formally exact discussion 



1^ in terms of the exact eigenstates in the presence of disorder. We will see that when 
the disorder potential in the two layers are similar, Vi{r) ~ Vr{r), the characteristics of 
the resonant tunneling change qualitatively when no magnetic is presented. The resonant- 
tunneling peaks are always sharp, no matter how strong the disorder; i.e., disorder no longer 
suppresses the divergence in the tunneling conductance. 

To leading order in — t — we may use the Fermi golden rule expression for the tunneling 
conductance |15[(-B = 0): 

■ ^(e/,r-e/3r)] (14) 

a/3 

where "^ai ^/3r are the exact eigenstates of the disordered 2DEG systems and the second 
equality holds to leading order in the voltage difference between the layers. When the 
electrons in different layers have identical disorder potentials Vi{r) = K-(r), the eigenstates 
are orthonormal {^ail'^iSr) = Sa,i3 and eai — €ar = Eoi — Ear, {E^i Eor are the energies at 
the subband edges). In this case Eq. (??) gives the tunneling conductance as: 




(15) 



A (5-function like sharp resonant tunneling peak occurs regardless of the amount of disorder 
presented, provided that Vi{r) = Vr{r). 
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Clearly, this exact result is not captured in the discussion of section V. In Feynman 
diagram language, we have to sum over diagrams which are important when the disorder 
potentials in the two layers are strongly correlated. When disorder is treated in the Born 
approximation the vertex-correction which yields a conserving approximation is shown in 
Fig.(l^); it is given by a sum of ladder diagrams with electron propagators of opposite 2DES's 
connected by the rungs of the ladder. The first diagram in Fig.( |^) is the 'bubble' diagram 
which gives the results of the previous section. 

Snhii^^) =^GkriiPn)Qkiiipn + iuJn) (16) 



Adding the vertex correction: 

^RJ^i'^^) = XI ^r{p, iPn)Gl{p, iPn + iu})Tf;{ipn, ipn + iuj) (17) 

where F satisfies the inhomogeneous integral equation, 

( iPn.iPn + i^) 

= 1 + V fr\{k - k')Qr{K', ipn)Qi{K\ ipn + iuj) 

K' 

■ Tj^,{ipn,ipn + i^) (18) 
Here the impurity-averaged disorder potential is: 



{Vrimir)) = ^ E /ir(p)e^^>-") (19) 
p 

To illustrate the physics , we consider the following short range model-potential calcula- 
tion: 

frrip) = Up) = 

Up) = xu' (20) 

where U^and A are constants. The case of identical impurity potentials in the two 2DES's 
corresponds to A = 1, and the uncorrelated layer limit in which the 'bubble' approximation 
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is appropriate corresponds to A = 0. This simple disorder-potential model allows us to solve 
the integral equation [Eq .(ll^)] very easily, since T is now independent of wavevector: 

T{ipn, ipn + ^t^) = 7775 (21) 

With the above expression for F, we can evaluate Eq. (|1^) by contour integration |p. After 
some standard manipulations we find that the conductance for c<j — and (3 ^ 00 (low 
temperatures) is : 

0+22 

G = — ^ E Re{Gf''{p, 0)Gr{p, 0)T,i-zO\ zO^) 

- G'f"(p,0)G'f"(p,0)^p^(^O^^Ot)} (22) 

Eq.(|2^) is easily evaluated numerically and the results for 5 = are shown in Fig.( ^ 
(The second term on the right hand side becomes negligible compared to the first for 
rep/h ^ 1. We see that, for A = 1, the resonant tunneling is singularly sharp even for 
a finite lifetime; this result is qualitatively different from the situation with A < 1 and 
demonstrates that the ladder captures the exact result discussed above. Numerical results 
with finite magnetic fields are shown in Fig.( ^); we see that while vertex corrections still 
reduce the effects of disorder the effect is quickly suppressed with increasing field. 

In typical double-layer samples the dopant layers are placed outside of each 2DES. In 
this case the disorder in each layer will come mainly from the closest dopant layer and we 
can expect the disorder potentials in the two layers to be weakly correlated. However in 
some circumstances it may be desirable to place the doping layer in the middle of the barrier 
separating 2DES's. In this case the disorder potentials in the two layers will be strongly 
correlated. The calculations in this section show that, even though doping in the center 
of the barrier is likely to decrease r substantially it will not effect tunneling between the 
2DES's at 5 = 0. m 
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V. SUMMARY 



In summary we have derived and evaluated expressions for the low temperature dc equi- 
librium tunneling conductance between parallel two-dimensional electron systems. Analytic 
results were obtained for the dependence of the tunneling conductance on the strength of 
an in-plane magnetic field and on the layer density difference between the 2DES's. The 
effect of uncorrelated disorder on the tunnehng conductance was discussed. Expression re- 
lating measures of the broadening of conductance curves by disorder to rep/h were given 
and the possibihty of using them to measure the electron lifetime in a 2DES was discussed. 
Comparisons of our theoretical results to the experimental results of Eisenstein et al. demon- 
strate that non-momentum-conserving tunneling events are remarkably infrequent in their 
samples. Finally we have shown that vertex corrections become important, especially at 
zero in-plane field, when strong correlations exist between the disorder potentials in the two 
2DES's. These vertex corrections are likely to be important in practice in samples where 
the doping layer is in the middle of the barrier separating the 2DES's. 
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FIGURES 

FIG. 1. Tunneling conductance vs. in-plane magnetic field strength at equal layer density, 

ni = Tlr = 1.5 X loHcm-2. The insert shows the displaced Fermi circles in the two layers. Momen- 
tum conserving tunneling events at the Fermi energy are allowed at the intersections of the two 
Fermi circles. The dashed line shows the result in the absence of disorder while the solid line was 
evaluated for a quantum lifetime r = Ibh/ep- Gq/A = 1.5 x 1.0'^{m*'^e'^t'^ /n'^h^ni). 

FIG. 2. Tunneling conductance vs. the ratio of electron densities in the two layers. 
{ni = 0.9 X lO^^cm"^ These results were calculated with rep/h = 15. The dashed line shows the 
result obtained in the absence of disorder at B = 0.4Tesla. Gq/A = 0.9 x l.(f{m*'^eH'^/Tr^h^ni). 

FIG. 3. Ladder diagrams for the vertex correction approximation for the current current re- 
sponse function. 

FIG. 4. Tunneling conductance versus layer density ratios at zero in-plane magnetic field for 
different degrees of correlation between the impurity potentials in the two 2DES's. As the cor- 
relation becomes strong (A 1) tunnelling occurs only at nearly equal densities even if there is 
substantial disorder in each layer. 

FIG. 5. Tunneling conductance versus layer density ratios at in-plane magnetic field B = O.AT 
for different degrees of correlation between the impurity potentials in the two 2DES's. The corre- 
lation does not change the tunnelling behavior qualitatively. 
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